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Abstract 



This talk serves as an introduction to higher twist effects in nuclei. 
£T) I We want to discuss how perturbative QCD can be applied to processes 

involving heavy nuclei by taking into account multiple scattering. 

1 Introduction 

en ■ 

Perturbative Quantum Chromodynamics (pQCD) has been established as an 
extremely successful theory to describe phenomena related to scattering reac- 
tions off hadrons at large momentum transfer. pQCD calculations combined 
with the available parameterizations of parton distributions enable us to ex- 
Ph! plain a large set of data including those from deep inelastic lepton nucleon 

scattering (DIS) I + N — ► I + X and the famous Drell Yan process (DY) 
N + N — >l+ + l~+X. 

pQCD is in fact a rather strict theory in the sense that no assumptions 
arising from any model enter. There are quantities of non-perturbative nature 
which cannot be described by perturbation theory, e.g. the bound states of 
QCD. Nevertheless there exists a quite rigorous way to separate perturbative 
(short range) and non-perturbative (long range) physics in a scattering reaction 
in a proper way. Let us consider DIS as an example. A factorization theorem |T| 
enables us to shift all non-perturbative physics into a set of well-defined, gauge- 
invariant (i.e. observable) and universal (i.e. process independent) quantities. 
These quantities can be expressed by matrix elements of parton operators 
between hadron states. 

In principle an infinite number of such matrix elements could enter our 
calculation and spoil the usability of the factorization theorem. But we can 
establish a hierarchy between the matrix elements in terms of an expansion 
in inverse powers of the momentum transfer. In more detail the quantity we 
expand in is X/Q where Q is the perturbative hard scale, e.g. the virtuality 



of the virtual photon in DIS, and A (for massless QCD) has to be some non- 
perturbative (hadronic) scale. The perturbative scale Q makes the coupling 
running and has to be large (at least a few GeV) in order to make perturbation 
theory applicable while A is of the order of Aqcd ~ 200 MeV. Under this 
circumstances X/Q <C 1 and this expansion, called the twist expansion, will 
do well. This is the definition of twist we will use here. 

The factorization theorem tells us that the left diagram in Fig. 1 gives 
the leading contribution to DIS in the twist expansion (and also in the a s 
expansion) .[] 



Fig. 1: Left: leading twist contribution to deep inelastic lepton nucleon scattering. The 
virtual photon emitted from the lepton (not shown) scatters off one quark inside the nucleon. 
Right: leading twist contribution to Drell Yan pair production in nucleon nucleon collisions. 
Each nucleon emits a quark or antiquark which annihilate into a virtual photon, which then 
decays into a lepton pair (not shown). 

Note that the leading twist contribution always consists of one hard scat- 
tering on the parton level, here the scattering of the photon off one quark 
from the nucleon. The non-perturbative part is described by a matrix element 
which encodes the process of taking one quark out of the nucleon and putting 
it back (in the complex conjugated graph). It is given by 



where q are quark field operators and P is the momentum of the nucleon. 
The Bjorken variable £ gives the momentum fraction of the nucleon that is 
carried by the quark in an infinite momentum frame. / 9 (£) is then just the 
famous quark distribution function in the nucleon. It has a parton model 
interpretation as the probability to find a quark with momentum fraction £ in 
the nucleon. This leading twist contribution to DIS is also called twist-2. The 
next order (next-to-leading twist) would involve matrix elements of twist-4, 

1 We are not going to discuss the subtleties arising from the fact that we have two expan- 
sions, one in X/Q and one in a s ~ In X/Q, which are obviously not independent. 
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i.e. 0(A 2 /Q 2 ).Q The cross section for DIS can then be written as a convolution 
of a parton distribution f a with a parton cross section ai +a — which itself is a 
series in a s — plus power corrections which arise from matrix elements which 
differ from parton distributions. 

ai +p = J2f a ®ai +a + o(jpj (2) 

For the Drell Yan process also a factorization theorem holds which gives the 
leading twist contribution as a simple annihilation of quarks and antiquarks 
to produce a virtual photon, see Fig. 1 (right). Note again that leading twist 
(twist-2) involves only one hard scattering on the parton level (the annihila- 
tion), while two matrix elements, one for each nucleon, are in the game. These 
matrix elements define exactly the same parton distributions as in DIS. This 
powerful property is called universality. Measuring a parton distribution in 
one process gives us predictive power for all other processes where this parton 
distribution enters. 



2 Multiple Scattering and Nuclear Enhancement 

What changes if we replace single hadrons by nuclei and look at e + A, p + A 
or A + A collisions? The factorization theorems should hold also here, but for 
many observables the picture of one hard scattering without rescattering on the 
parton level seems not to be the dominant one. E.g. for Drell Yan we expect 
initial state interactions when the quark and the antiquark participating in the 
annihilation have to traverse a large piece of nuclear matter. In the framework 
of pQCD multiple scatterings are exactly higher twist corrections to the single 
scattering (leading twist) process. 

At a given scale Q 2 the power corrections by higher twist effects are en- 
hanced because of the nuclear size. This enhancement cannot stem from the 
hard (short range) part of the cross section, but must arise from the non- 
perturbative (long range) matrix elements. We will demonstrate later how the 
matrix elements can test the nuclear size. 

It was an idea mainly advocated by Luo, Qiu and Sterman 0] that, e.g. on 
the twist-4 level, there are matrix elements which scale like \ 2 A l l 3 /Q 2 where A 
is the mass number of the nucleus, whereas there are others which are not sen- 
sitive to the nuclear size and behave just like X 2 /Q 2 as usual. For large nuclei 
we can assume that A 1 / 3 ^> 1 and conclude that only those higher twist cor- 
rections, which show an additional scaling with the nuclear size, called nuclear 

2 Twist-3 vanishes for unpolarized scattering. 



enhanced, are important and we safely omit all others. We rearrange our twist 
series in such a way, that we have an expansion in powers of A 2 A 1 /' 3 /Q 2 , i.e. 
for each power of Q~ 2 we have an additional power of A 1 ^ 3 (maximal nuclear 
enhancement). Terms with less powers in A 1 ^ 3 than in Q~ 2 are subleading. 

Twist corrections normally are a tricky business and very difficult to handle. 
Starting from the parton distribution in Eq. (1) we can e.g. generate higher 
twist matrix elements just by inserting covariant derivatives acting on the 
parton fields. Unfortunately the number of matrix elements contributing at a 
certain level of twist can be large and in addition most of these matrix elements 
have no probabilistic interpretation in a parton picture. 

This changes dramatically if we demand that a matrix element should be 
sensible to the size of an extended medium. Extended here means that the 
extension L should be much larger than the confinement radius Rq, L ^> R Q , 
what is true for large nuclei. One can show that only matrix elements with 
additional pairs of parton operators can contribute. For twist-4 schematically 
we have e.g. (P\ (qq)(FF) \P) or (P\ (FF)(FF) \P), where color, spinor and 
vector indices for each pair are contracted in the same way as in the parton 
distribution respectively. F here is a gluon field strength. Therefore these ma- 
trix elements are limited in number and have a straight forward interpretation 
as correlators between partons. 

In Fig. 2 we give an example for single and double scattering in the DY 
process at large transverse momentum. Here the quark (antiquark) from the 
single hadron has the opportunity to scatter off an additional gluon from the 
nucleus before annihilating with an antiquark (quark). 
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Fig. 2: Possible diagrams for single scattering (left) and double scattering (right) contribut- 
ing to DY pair production in p + A collisions at high transverse momentum. Note that both 
diagrams are next-to-leading order in a s respectively, since we radiate an extra gluon in 
order to generate finite transverse momentum. 

The onset of nuclear enhancement can be seen in p + A experiments. The 
best known example is the famous Cronin effect || |j], which is a nuclear en- 
hancement of the production cross section of Drell Yan pairs, J/ips or other 
particles at large transverse momentum, whereas at low transverse momen- 
tum enhancement id missing. The enhancement is a clear sign of multiple 



scattering. Another sign is the transverse momentum broadening of jets or 
Drell Yan pairs. Fermilab data indicate a rise of the ratio of second and first 
moment of the transverse momentum spectrum which is compatible with an 
A 1 ^ behaviour [|J. 

Let us emphasize that single scattering contributions in principle scale with 
A, i.e. just they scale with the volume the reaction can take place. The nuclear 
enhanced double scattering picks up an additional power of the nuclear radius 
and therefore scales like v4 4//3 and so on. With this mind one would expect 
an overall enhancement of all cross sections but experimentally this is clearly 
not the case. Indeed most cross sections are not enhanced like e.g. DY at 
low transverse momenta. If pQCD intends to explain the enhancement by the 
introduction of new matrix elements, it should also account for the absence 
of enhancement in other observables. Again this works fine. We will see that 
interference seems to be very important. This can spoil the enhancement in 
certain observables. 

3 An example: double scattering in Drell Yan 

Let us discuss the Drell Yan process in p + A as an example. We would like to 
calculate the effect of double scattering at large transverse momentum || ^, § . 
The DY cross section with full kinematic dependence on Mass Q, transverse 
momentum q± and rapidity y of the lepton pair and on the angular distribution 
of the lepton pair can be decomposed into four structures parametrized by four 
helicity amplitudes Wtl, Wl, and Waa- 

_^_(w T1 (l + cos^) + ^(i-5co^) + 

W A (sin 26 cos 0) + W AA (sin 2 9 cos 20) ) (3) 

If we integrate over the angles, only the term proportional to Wtl contributes 
to the cross section. Fig. 3 shows two typical examples for graphs for twist-4. 
At this order of twist we not only have to take into account double scattering 
but also interference of single and triple scattering (which we will include into 
the term double scattering in the following since the matrix elements are the 
same) . 
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Fig. 3: Two diagrams contributing to twist-4 at large transverse momentum. The left 
diagram shows an ordinary double scattering with a quark gluon pair from the nucleus on 
the left and right (complex conjugated) part of the diagram. The right diagram shows a 
possible interference between triple scattering (two gluons and one quark on the left part) 
and single scattering (one quark on the right side). 

It turns out that there are two different contributions to double scattering. 
Both partons from the nucleus can be hard (double hard scattering) or one can 
be hard and one soft (soft hard scattering). The final formulae for the helicity 
amplitudes Wi including the partons a and b from the nucleus and parton c 
from the single nucleon are 

W t DH ' ab+C = f ^rM&T™^, x h )Hr' ab+c (x a , x h , 6), (4) 



w SH, ab+ c = 

JB 42 

for double hard and soft hard scattering respectively. f c are the usual parton 
distributions for the nucleon, H ab+C is the cross section on the parton level 
(the hard part of the cross section) and the are the new matrix elements 
describing two partons in the nucleus. 

Let us shortly discuss where the distinction between soft hard and double 
hard arises from. We have to integrate over longitudinal momenta of parton 
lines connecting soft and hard parts of the diagram and we have propagators in 
the hard part which can provide poles in the integrands. An example is given 
in Fig. 4. Both propagators marked by circles give poles for the integration of 
the gluon momentum fraction x. 




Fig. 4: Soft pole (left) and hard pole (right) indicated by circles. 

The upper pole fixes x to a large (finite) value Xhard, whereas the lower pole 
gives x so f t ~ Kj_/S ~ where K± is the intrinsic transverse momentum of the 
gluon. Carrying out the pole integration gives a sum of both residues: 

M ~ / dx ( 1 -—t-) ( F(x, k t , 6) 

J \x- x soft + \x- x hard + iej 

\ %soft -^hard / \ -^soft %hard 

— M so f t ^ard — ^haTd,hard 

Note the minus sign between both terms. To get the cross section we have to 
square the term above. The hierarchy Xhard 3> Xsoft can only be ensured as long 
as the transverse momentum q± is large. In this case the interference between 
both residues is negligible and the cross section is given just by the sum of 
squares Mg f thard 

+ M Lrd hard- The first term is the soft hard contribution 
(soft gluon) the second term is the double hard contribution (hard gluon). 

However if q\ -C Q 2 then Xhard starts to approach x so f t and both residues 
start to cancel each other. This is the explanation of the important state- 
ment already given above: for small transverse momenta the double scattering 
contribution undergoes destructive interference and there is no nuclear en- 
hancement from double scattering in that kinematic region. 

Let us now discuss the twist-4 matrix elements. The double hard matrix 
elements depend on the momentum fractions £ and x of both hard partons. 
The quark gluon correlator e.g. is 

tS?&x) = ^/d^^e(, r -%)e(-c) (r) 

^PtH^PtiH~^)\ F" + {zl)F + u (z^)q(0)-y + q(^) \Pi) ■ 
The same correlator for soft hard scattering reads 

T q f(0 = J dz^^e(zT - z^)e(- zl ) (8) 




e ^V_ ( Pl | F"+{zl)F+ ^miMzi) \Pi) ■ 



It depends only on the momentum of the hard parton. The 0-functions en- 
sure causality. PQCD is not able to predict these matrix elements from first 
principles and they are not measured up to now. In order to make any numer- 
ical statements about the size of the nuclear enhanced corrections we have to 
rely on models at this stage. However the important point is, we can extract 
the scaling with the nuclear size. To do this we have to take into account 
the colour structure of the operators, the oscillating exponential factors and 
the boundaries of the integrals. An analysis gives that both matrix elements 
above have one free integration which can test the extension of the nucleus. 
Pictorially this is the distance between the two parton pairs qq <-> FF. From 
that we assume that the matrix elements can be modelled by ordinary parton 
distributions by setting 



C and A 2 are normalization constants of the order of Aqcd- 

At the moment our goal must be to gather experimental information about 
the new matrix elements and therefore our study for RHIC |7|, || intends to look 
for observables where this can be done. Independent of the models which one 
plugs in for the matrix elements one can establish some results: Double hard 
scattering has trivial angular dependence like the leading twist, leading a 5 
process. It seems to confirm the picture of two independent binary collisions, 
first q+g — > q+g and then q+q — > 7*. Confer the right diagram of Fig. 4: the 
pole splits the diagram into two independent subprocesses. Also the so called 
Lam Tung sum rule 2Waa = Wl, a long standing leading twist prediction 
[P], is respected by double hard scattering. On the other hand the results on 
soft hard scattering are more complicated and violate the Lam Tung relation. 
Another important point is that the leading twist result for is nearly zero 
(it vanishes for p + p from symmetry reasons), while there is a non vanishing 
twist-4 contribution. These are the most important model independent results 
which could allow a glance at the new matrix elements in forthcoming data. 
At the end of this section we show two results for proton gold collisions at 
RHIC. In Fig. 5 the rapidity distributions of the amplitudes Wtl an d Wa are 
given. Note that at RHIC energies soft hard seems to be dominated by the 
double hard process. 
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Fig. 5: Results for 100 GeV Au + 250 GeV p collisions at RHIC. Cross sections for helic- 
ity channels Wtl (left) and Wa (right) as a function of rapidity y at Q = 5 GeV and 
q± = 4 GeV. Single scattering (solid lines), double hard (dashed) and soft hard (dot 
dashed) for two different choices of parton distribution parameterizations (CTEQ3M and 
CTEQ3M+EKS98) are shown. 



4 Further topics 



The transverse momentum broadening of DY pairs was also explained recently 
by X. Guo within this formalism |HJ. For this observable one is sensible to 
small transverse momenta and one has to calculate graphs like in Fig. 6 (left) 
with no additional parton in the final state. The ratio of second and first 
moment of the transverse momentum spectrum is surprisingly easy 

(ql) 4ir 2 a s f c 
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(ql) 3 f c ®f a 3 

with partons a and c being quark and antiquark or vice versa. The last equal 
sign holds for our model for the soft hard matrix elements. One may note that 
the pQCD result indeed rises with A 1 / 3 compatible with the data. 

In the low q± regime one would like to sum the effect of multiple scattering 
which would be a resummation of higher twist. Diagrams like in Fig. 6 (right) 
introduce matrix elements which are correlators with an arbitrary number n 
of gluon fields (P\(qq)(FF) n \P) . If we extrapolate the model above we can 
establish connections between matrix elements of different n and are able to 
carry out the summation. 





Fig. 6: DY graphs with scattering off one (left) and n (right) additional gluons from the 
nucleus. 



For the transverse momentum spectrum we obtain a shift and smearing. 
The shift is given by Aii 2 a s A l ^\ 2 /3. Some other topics addressed recently 



include the medium modification of fragmentation functions [11] which can 
explain energy loss and the medium modification of parton distributions via 
higher twist modifications internal to the nucleus [0 which can account for 
shadowing. 

Let us close with some remarks about the limitations of pQCD fl2fl . It 
is well known that QCD factorization breaks down in Drell Yan beyond the 
order of twist-4 Practically this means that there are non-factorizable 

contributions that make the twist expansion no longer reasonable beyond twist- 
4. However the nuclear enhancement argument saves us for p + A collisions 
since it allows us to stay at twist-2 for the single proton. In that case we can 
go to arbitrary twist for the nucleus, like we can do also for e + A. For A + A 
however the problem cannot be waived and we are limited to twist-4 accuracy, 
which this still gives us the leading medium effect. 
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